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INTRODUCTION 


The  requirement  to  predict  satellite  trajectories  accurately 
depends  very  critically  on  detailed  knowledge  of  upper  atmospheric 
density.  The  atmospheric  model  currently  being  used  at  the  Satellite 
Control  Facility  (SCF)  in  Sunnyvale  has  errors  in  it  which  can  lead  to 
corresponding  position  errors  which  are  unacceptably  large  at  the 
end  of  a one -day  period.  In  addition,  the  long-term  lifetime  predic- 
tions based  on  this  model  have  little  validity.  The  situation  is  most 
critical  for  high  inclination  orbits  at  low  altitudes.  The  errors  are 
expected  to  become  more  severe  as  solar  maximum  is  approached 
around  1980. 


The  model  currently  in  use  at  SCF  is  called'DENSEL.  This 
model  is  essentially  a juxtaposition  of  the  ARDC  1959  (Air  Research 
and  Defense  Command)  model  (Mizner  et  al.  , 1959)  and  the  first  thermo- 
spheric density  model  produced  by  Jacchia  (I960).  The  simplicity  of  the 
'DENSEL  model  makes  it  useful  for  production-type  calculations  such 
as  orbit  determinations.  However,  since  it  doesn't  include  any  of  the 
known  atmospheric  variations  other  than  the  diurnal  and  solar  cycle 
effects,  its  accuracy  is  rather  poor.  In  addition,  recent  studies  by 
Sharp  et  al.  (1978)  have  shown  that  the  diurnal  variation  in  'DENSEL  and 
other  models  derived  from  satellite  drag  data  is  incorrectly  modeled, 
especially  at  low  altitudes.  Also,  since  there  is  no  provision  for 
input  to  the  model,  the  density  predictions  of  'DENSEL  are  too  low 
during  a magnetic  storm.  As  a consequence  of  these  errors  in  'DENSEL, 
the  right  combination  of  K^,  month  of  the  year,  and  local  time  could  lead 
to  a situation  where  'DENSEL  would  predict  a density  that  was  too  low  by 
more  than  a factor  of  2. 
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DISCUSSION 


In  the  last  18  years  enough  progress  has  been  made  in  thermo- 
spheric physics  and  upper  atmospheric  modeling  to  justify  replacing 
'DENSEL  with  a more  accurate  model.  In  order  to  select  a replace- 
ment model  for  'DENSEL,  extensive  tests  were  run  on  five  current 
upper  atmospheric  density  models  and  'DENSEL  using  the  neutral 
density  data  base  compiled  at  the  Space  Sciences  Laboratory  of  The 
Aerospace  Corporation.  This  data  base  currently  consists  of  over 
40,  000  measurements  of  neutral  atmospheric  density  at  5 km  inter- 
vals in  the  range  130-400  km.  The  data  were  obtained  from  three 
different  types  of  instruments:  cold  cathode  ion  gauges,  capacitance 
manometers  and  calorimeters.  The  former  two  instruments  are 
described  by  Rice  et  al.  (1973).  The  data  reduction  procedures  for 
constructing  the  data  base  have  been  described  by  Rice  and  Sharp 
(1977). 

The  five  other  models  tested  were  MDAC  (Olson  et  al.  , 1975) 
MSIS  (Hedin  et  al.  , 1977a,  b),  J-65  (Jacchia,  1965),  J-71  (Jacchia, 

1971)  and  J-77  (Jacchia  1977).  A comprehensive  description  of  all 
these  models  with  the  relative  advantages  and  limitations  of  each 
has  been  provided  by  Hickman  et  al.  , (1978).  Two  measures  of  a 
model's  effectiveness  were  employed  in  this  evalxxation,  the  average  per- 
centage error  and  the  RMS  error.  The  average  absolute  precentage 
error  is  simply: 
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Where: 


p = Measured  Density 

nri  = Density  Predicted  by  the  Model 

N = Number  of  Points 


The  RMS  error  is  defined  as: 


C 


In  the  RMS  error  calculation  any  offset  between  the  model  and  the  data 
base  is  removed  and  the  larger  amplitude  errors  are  weighted  more 
heavily. 

In  Table  1 is  presented  the  results  of  this  test  divided  into  Z 
altitude  regimes,  130-195  and  200-400  km.  Since  MSIS  performed 
the  best  in  the  low  altitude  regime  where  density  predictions  are  most 
critical,it  was  selected  as  the  model  to  replace  'DENSEL. 

The  MSIS  (Mass  Spectrometer  and  Incoherent  Scatter)  model 

makes  use  of  measurements  of  upper -atmospheric  compositions  from 

Mass  Spectrometers  on  five  satellites  (AE-B,  OGO-6,  San  Marco-3, 

Aeros-A  and  AE-C)  and  neutral  temperatures  inferred  from  incoherent 

scatter  measurements  at  four  ground  stations  (Arecibo,  Jicamarca, 

Millstone  Hill,  and  St.  Santin).  The  data  set  covers  the  time  period  from 

1965  to  1975,  essentially  a complete  11  year  solar  cycle.  There  are  two 

liJTiitations  in  MSIS  which  should  be  corrected  to  enhance  the  accuracy  of 

the  model  predictions.  The  first  is  the  fact  that  MSIS  has  no  K input.  It  models 

P 

the  corpuscular  heat  source  simply  by  using  the  daily  Ap  index.  Since  the 
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TABLE  1 


All  Data 

Number  of  Points  40,  261 

130-195  km 

25, 639 

200-400km 

14. 622 

Average  Absolute  Percentage  Error 

'DENSEL 

22.  3 

28.  9 

26.2 

MSIS 

13.  5 

12.  6 

15.  8 

J-71 

14.  2 

13.  5 

16.  3 

J-77 

15.  4 

15.  1 

17.0 

J-65 

18.  0 

14.  7 

30.4 

MDAC 

28.  8 

23.  8 

30.  6 

Average  RMS  Percentag 

e Error 

'DENSEL 

23.4 

17.  0 

28.  7 

MSIS 

16.9 

14.  2 

18.  2 

J-71 

16.  7 

15.  7 

18.  8 

J-77 

17.  8 

17.  5 

18.  8 

J-65 

18.  3 

16.  1 

18.9 

MDAC 

21.0 

19.  1 

21. 8 

thermosphere  is  known  to  respond  to  changes  in  with  a delay  time 

of  about  6 hours;  a correction  algorithm  has  been  added  in  a form 
which  best  fits  The  Aerospace  Data  Base.  For  ^ 4.  3 the  correction 
is  given  by 


D (6)  = D (6)  X r ( K +1.0  1 

- p 


where  D (6)  is  the  output  density  of  the  basic  model.  For  K > 4.  3 

P 

the  correction  is  given  by 


r f K .0.4 

D ,6)  . D (6)  , [ ( .5^  ) ] 


For  Kp  = 4.0  the  correction  factor  from  the  first  equation  is  1.083 
and  for  = 4.  3 the  correction  factor  from  second  equation  is  1.085; 

for  Kp  = 8 the  correction  factor  is  1.391.  The  application  of  this 
correction  algorithm  yields  a model  known  as  MSIS  (MOD  1 ). 

The  second  correction  is  applicable  to  the  altitude  dependence 
of  the  model  output.  MSIS  in  general  predicts  densities  which  are  too 
low  at  low  altitudes  and  too  high  at  high  altitudes  when  compared  with 
The  Aerospace  Data  Base  and  the  Jacchia  Models.  The  following  correc' 
tion  algorithm  compensates  for  this  effect 


D (6)  = D (6)  X r-98  + ( 


185-H 

1200 


) ] 


where  H is  the  altitude  in  km. 
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At  130  km  the  correction  factor  is  1. 025  and  at  300  km  the  correction 


factor  is  .884.  The  addition  of  this  correction  algorithm  to  MSIS 
(MOD  1)  yields  a model  known  as  MSlS  (MOD  2).  The  accuracy  analysis 
for  these  two  modifications  is  presented  in  Table  2 along  with  the 
original  results  for  MSIS  and  'DENSEL,  repeated  for  convenience. 

As  can  be  seen  from  Table  2,  the  accuracy  difference  between 
'DENSEL  and  MSIS  (MOD  2)  is  significant  but  not  overwhelming,  so 
an  expensive  systems  software  change  was  deemed  inappropriate.  Con- 
sequently, a set  of  correction  tables  was  generated  which  provide  the 
ratio  of  MSIS  (MOD  2)  to  'DENSEL  for  a number  of  input  parameters 
designed  for  wide  coverage,  good  resultion  and  ease  of  use.  The 
computer  Fortran  Code  is  contained  in  Appendix  A and  a portion  of  the 
printout  is  contained  in  Appendix  B.  Copies  of  subroutines  'DENSEL 
and  MSIS  are  available  upon  request. 

The  format  of  the  computer  printout  is  set  up  in  such  a way  so 
the  user  can  rapidly  enter  the  appropriate  input  parameters  and  find  the 
necessary  density  correction.  The  first  column  gives  the  month  of  the 
year.  All  12  months  are  provided  since  'DENSEL  has  no  semiannual  vari- 
ation. The  day  of  the  month  is  assumed  to  be  the  15th.  The  next  column 

is  the  altitude  in  nautical  miles.  The  tables  have  been  provided  in  two 

i 

sets,  one  from  65  to  95  NM  and  the  second  set  from  100  to  150  NM,  both  in 
10  NM  steps.  The  third  column  is  latitude  with  four  values  provided 
(0°,  20°,  40°,  60°N).  The  fourth  column  is  local  time  with  six  values 
provided  (2.0,  6.0,  10.0,  14.0,  18.0,  22.0).  The  fifth  column  is  the 
■ 1 0.  7 cm  flux  for  the  day  prior  to  the  day  the  prediction  is  to  be  made 

with  four  values  provided  (80,  120,  160,  200).  The  sixth  column  is  the 
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TABLE  2 


Average  Absolute  Percentage  Error 

•DENSEL 

22.3 

MSIS 

13.5 

MSIS  (MOD  1 ) 

13.1 

MSIS  (MOD  2) 

12.1 

Average  RMS  Percentage 

Error 

'DENSEL 

23.4 

MSIS 

16.9 

MSIS  (MOD  1) 

16.4 

MSIS  (MOD  2) 

15.4 
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Kp  value  which  occurred  six  hours  prior  to  the  time  for  which  the  pre- 
diction is  to  be  made  with  five  values  provided  (0,  1,  A,  U,  8),  The 
eighth  column  is  the  ratio  of  the  MSIS  (MOD  2)  density  prediction  to  the 
'DENSEL  prediction.  The  next  column  is  labeled  B-Factor  (which  is 
short  for  "B-Factor  Correction")  since  this  ratio  should  be  proportional 
to  the  B-Factor  which  is  currently  being  generated  from  observation  of 
orbit  precession  and  knowledge  of  offsets  and  errors  in  'DENSEL.  The 

last  two  columns  give  the  actual  densities  for  both  MSIS  (MOD  2)  and 

3 

'DENSEL  in  grams  per  cm  for  reference.  For  best  possible  accuracy 
in  using  these  tables,  some  interpolation  will  be  necessary. 


CONCLUSION 


Two  sets  of  correction  tables  have  been  provided  to  the  Satellite 
Control  Facility  in  Sunnyvale,  California  which  will  aid  in  satellite 
tracking  and  navigation.  These  tables,  one  for  low  altitudes  and  the 
other  for  high  altitudes,  incorporate  a modified  version  of  the  MSIS 
Model  and  provide  the  user  with  a ratio  of  MSIS  (MOD  2,)  to  'DENSEL 
for  a variety  of  input  parameters.  These  tables  shall  eliminate  the 
need  to  make  ad  hoc  corrections  to  the  'DENSEL-derived  B-Factors 
since  these  systematic  corrections  are  derived  from  the  best  currently 
available  atmospheric  model.  Thus,  short-term  ephemeris  predictions 
and  long-term  lifetime  estimates  should  be  considerably  more  accurate 
than  they  are  at  present. 
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APPENDIX  B (Continued) 


THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver> 
satility  and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel> 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans> 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory;  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  iiTi^Wted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  ^pli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  piienomena,  including  plasma  electromagnetics;  quant^im  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting. superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials;  metal 
matrt3rTompo8ite8~an9~new~formFM  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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